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1 INTRODUCTION 



A number of independent cosmol ogical observations, ranging from 
type- la supernovae (SNIa) (e.g. iRiess et al.ll2007l : |Perlmutter et al.l 
19991). to Cosm ic Mic rowave Background (CMB) anisotropies (e.g. 
Koma tsuet alj|2oTll; lLarson et alj|201l | ) and Large S cale Struc- 
ture (LSS) (e.g. IVikhlinin et alj|2009l : iReid et al]|20ld) . convinc- 
ingly show that the Universe is undergoing a phase of accel- 
erated expansion. One of the main challenges of modern cos- 
mology is in fact understanding the source of such acceler- 
ation. To this purpose, a number of models have been pro- 
posed that modify the two pillars of modern cosmology, gen- 



ABSTRACT 

We present forecasts on the capability of future wide-area high-sensitivity X-ray surveys 
of galaxy clusters to yield constraints on the parameters defining the Dark Energy (DE) equa- 
tion of state (EoS). Our analysis is carried out for future X-ray surveys which have enough 
sensitivity to provide accurate measurements of X-ray mass proxies and Fe-line based red- 
shifts for about 2 x 10 4 clusters. We base our analysis on the Fisher Matrix formalism, by 
combining information on the cluster number counts and power spectrum, also including, for 
the first time in the analysis of the large scale cluster distribution, the effect of linear redshift 
space distortions. This study is performed with the main purpose of dissecting the cosmologi- 
cal information provided by geometrical and growth tests, which are both included in the anal- 
ysis of number counts and clustering of galaxy clusters. We compare cosmological constraints 
obtained by assuming different levels of prior knowledge of the parameters which define the 
relation between cluster mass and X-ray observables. This comparison further demonstrates 
the fundamental importance of having a well calibrated observable-mass relation and, most 
importantly, its redshift evolution. Such a calibration can be achieved only by having at least 
~ 10 3 net photon counts for each cluster included in the survey, with sufficient angular resolu- 
tion. We show that redshift space distortions in the power spectrum analysis carry important 
cosmological information also when traced with galaxy clusters. We find that the DE FoM 
increases by a factor of 8 when including the effect of such distortions. Besides confirming 
the potential that large cluster surveys have in constraining the nature of Dark Energy, our 
analysis emphasizes that a full exploitation of the cosmological information carried by such 
surveys requires not only a large statistic but also a robust measurement of the mass proxies 
and redshift for a significant fraction of the cluster sample, which ought to be derived from 
the same X-ray survey data. This will be possible with future X-ray surveys, such as those en- 
visioned with the Wide Field X-ray Telescope, with an adequate combination of survey area, 
sensitivity and angular resolution. 

eral relat ivity and the standard m odel of fundamental interac- 
tions (e.g. lSilvestri & TrodderfeOOl and references therein). Mod- 
els that modify the latter and, therefore, the energy-momentum 
tensor in the Einstein equation are , for example, the scalar 
field models like quintessence (e.g. ICaldwell & Kamionkowski 
[2009|;lDoran & Wetterichl2003l;lziatev et alJl999l:lRatra & Peebles 



1988). k-essence (e.g. iTsuiikawaj |2010|: Mukohyama & Randall 
|2004) . c oupled Dark Energ y (e.g. lAme ndola 2000) and Chapligin 
gas (e.g. lBento et alj |2002). Models that modify general relativity 
can produce the cosmic acceleration without including a Dark En- 
ergy (DE) component, but they should also satisfy stringent con- 
straint from local (e.g. solar system) tests of gravity. Examples of 
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such models are the branewo rld models, like DGP (e.g. Dvali et al. 
200fj| ; )Movahed et alj 120091). f(R) theories (ej;. IHu & Sawickl 
2001 

2010T) and scalar-tensor theories (e 



Brax et al.l 120081 : lApplebv et alj l20ld; ISotiriou & Faraonil 
g. lSkordisll2009h . 
The first model proposed to explain cosmic acceleration was 
based on the introduction of a cosmological constant A that can 
be thought as a fluid with negative pressure and equation of state 
(EoS) which is constant in space and time, p = wpc 2 with w = -1. 
However, a cosmological constant able to drive the accelerated ex- 
pansion leads to the well known problem that A needs to be so tiny 
with respect to any natural energy scale, that there is no theoret- 
ical justification for it. Thus, a plethora of models, characterized 
by different p arametrizations o f the DE EoS evolution, have been 
proposed (e.g. IWetterichl2004l) . In principle, it is possible to distin- 
guish among different DE models by combining different cosmo- 
logical probes, both based on the geometry of the universal back- 
groun d and the growth of density perturbations (e.g. jAlbrecht et all 
l2006h . 

In this context, clusters of galaxies have long been recognized 
as potentially powerful pro bes of the nature of DE and cosmologi- 
cal models in general (e.g. lAllen et al.l201 lliLombriser et alj2010l : 
iManera & M ota 2006, and references therein). Clusters of galaxies 
provide cosmological information in a number of different ways. 
The evolution of the cluster space density depends on cosmologi- 
cal parameters through both the linear growth rate of density per- 
turbations and the redshift dependence of the volume element. The 
large-scale clustering of galaxy clusters is also sensitive to cosmo- 
logical parameters, through the growth rate of perturbations, which 
affects both the bias parameter and the redshift-space distortions, as 
well as by sampling the shape of the underlying Dark Matter (DM) 
power spectrum over a broad range of wavenumbers. As of today, 
relatively small samples of ~ 100 X-ray selected clusters, originally 
identified out to z — 0.8 by the ROSAT satellite and then followed- 
up by Chandra to obtain robus t mass estimates, have provided in- 
terest ing DE constraints (e.g.. IVikhlinin et "all 120091 : iMantz et al.1 
l2010l) . which complement and agree with tho se from CMB and 
SNIa observations (see also lRapetti et al.|[2005T) . More recently, by 
using a complete sample of < 1000 nearby cl usters at z~ 0.2, iden- 
tified in the ROSAT All-Sky Survey (RASS ; iTruempe j| 19931) . the 
large-scale powe r spectrum has been constructed o ver a fairly wide 
scale range (e.g., Bala guera-Antolmez et alj201 lh . 

While the first Sunyaev-Zeldovich (SZ) surveys have 
now started producing cluster sa m ples of similar sizes, (e.g.. 
Marriage & the ACT Team! 1201 ll; ("Williamson & the SPT Team! 
2011; Planck Collaboration 2011), the next generation of X-ray 
(e.g., eROSIT41 WFXT0) and optical (e.g., DES^ EUCLIE0) 
surveys are expected to increase by orders of magnitude the 
number of galaxy clusters, further extending the redshift range 
over which they trace the growth of cosmic structures. Such large 
cluster surveys have the potential of placing very tight constraints 
on different classes of DE models, possibly finding signatures of 
departures from the standard ACDM predictions. Several studies 
have dealt with constraints of DE models from future cluster 
surveys focusin g on the impact of uncertainties in cluster mass 



estimates (e .g. iBattve & Welled 120031; iMaiumdar & Mohr 
iLima & H j 12001 120071 : ICunhal |2009| ; ICunha & Evrard 



1 http://www.mpe.mpg.de/heg/www/Projects/erosita/index.php 

2 http://www.wfxt.eu/ 

3 http://www.darkenergysurvey.org/ 

4 http://sci.esa.int/science-e/www/area/index.cfm?fareaid=102 



2004; 



2010 



iBasilakos etafl koiO). All these analyses generally assume that, 
when a cluster is detected and included in a survey, the observable 
(i.e., X-ray luminosity, optical richness) on which the detection is 
based can be related to the actual cluster mass through a suitable 
relation, whose functional form is assumed to be known and 
depends on a set of additional parameters. 

A more conservative approach would instead require that for 
all clusters included in a survey detailed follow-up observations are 
carried out to calibrate suitable and robust mass proxies. Examples 
of suc h mass proxies in X-ray surveys include the total gas mass 
(e.g. IMantz et alJl2009h or the p roduct of gas mass a nd tempera- 
ture, the so-called Yx parameter (Kra vtsov et alj|2006l) . which can 
be computed when a relatively large number of photons is avail- 
able. Clearly, while measuring flux for an X-ray extended source 
requires only ~ 50 photons, or less for missions with low back- 
ground, the measurement of robust mass proxies requires at least 
~ 10 3 photons. 

In this paper, we will derive forecasts for constraints on the 
EoS of DE models from future X-ray surveys. We envision that 
these surveys will be carried out with a telescope with high-enough 
sensitivity to readily provide robust measurements of mass proxies 
and Fe-line based redshifts for ^ 10 4 clusters, which are all char- 
acterized by "C/iandra-quality" data, thus avoiding the need of ex- 
ternal and time-consuming follow-up observations. For example, 
the Wide Field X-ray Telescope (WFXT) concept was developed 
to meet such a requirement. This telescope combines large collect- 
ing area with a large field-of-view and sharp point spread function 
(PSF) approximately constant over the entire field of view (e.g., 
Giacconi et al] 120091 ; iMurrav & WFXT Team! l20ld; iRosati et al.l 
201 ll) . 



We adopt the specifications of the WFXT surveys, as an ex- 
ample of next generation X-ray cluster surveys, and we compute 
cosmologica l forecasts using the well-established Fisher Matrix ap- 
proach (e.g.. |Pod elson 2003]), to combine information from clus- 
ter number counts and large-scale clustering. We will quantify the 
constraints expected on DE models and their dependence on the 
knowledge of the relation between X-ray observables and cluster 
mass, for a range of survey strategies (i.e. depth vs. sky coverage). 
In the course of this study, we will discuss how number counts and 
the power spectrum of the large-scale distribution of clusters con- 
vey cosmological information. 

As a novel contribution in this paper, we will show how the 
detection of Baryonic Acoustic Oscillations (BAOs) and redshift 
space distortions (RSDs) on cluster scales can significantly con- 
tribute to constrain cosmological parameters, similarly to a number 
of previo us studies based on the large scale distribution of galax- 
ies (e.g.. iGuzzo et al.1 [200sl ; iRassat et alj 120081 ; IStril et al.ll201ol; 
IWang et alj|2010h . 

This paper is structured as follows. In Section 2, we briefly 
describe the parametrizations of the DE EoS that we use, and de- 
scribe our Fisher-Matrix approach for cluster number counts and 
power spectrum. In Section 3, we first describe the characteristics 
of the cluster survey, and derive the forecasts on the constraints 
on DE EoS parameters. In this section, we will also quantify the 
impact that uncertainties in the scaling relation between X-ray ob- 
servables and cluster mass have on such constraints. We discuss our 
results and present our conclusions in Section 4. 
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2 CLUSTERS AS DARK ENERGY PROBES 



2.1 Cluster number counts 



The relevant parameters of our analysis are the power spectrum nor- 
malization, cr 8 , the matter density parameter, Q,„, and the specific 
parameters defining the DE EoS. 

The reference analysis is carried out for the s tandard 
param etrization of the DE EoS, originally proposed by iLinderi 
(l2003h . 



The information Fisher Matrix (FM) is defined as 



w(a) = wo + w a (l — a) , 



(1) 



where a is the cosm ic expansion factor (see also 
IChevallier & Polarskil l200ll) . This parametrizati on has been 
used i n the D ark Energy Task Force reports (DETF; Albrec hFet al.l 
l200dl2009l) to assess the constraining power of different cosmo- 
logi cal experiments . 

lAlbrecht et alJfeOOeh presented forecasts on the constraints on 
the Wo a nd w a parameters fr om redshift number counts of cluster 
surveys. iMantz et alj d2010h derived constraints on these param- 
eters from the observed evolution of the cluster X-ray luminos- 
ity function , using a combination of nearby clust e rs selected from 
the R ASS jTruemped [l993l : lEbeling et all 1 19981 ; Bbhringer et al. 
|2004) and medium- distant clusters select ed from the (RASS based) 
MACS survey (e.g. lEbeling et al.ll2010f) . An update on the con- 
str aints available at pre sent on these parameters has been presented 
bv lKomatsu et all J201 lh using a combination of the 7-year WMAP 
CMB data, SN-Ia, Big-Bang Nucleosynthesis results and BAOs 
traced by the large-scale galaxy distribution. 

Furthermore, we will also assess the constraining power of 
X-ray cluster surveys for the class of quintess ence models, called 
Early Dark Energy (EDE) jWetterichl r2004h . In these models, 
DE drives not only the accelerated expansion of the Universe at 
relatively low redshift, but also provides a non-negligible con- 
tribution at early times, i.e. before the last scattering surface 
jPoran & Wettericrj|2003h . A parametrization of a class of EDE 
models has been proposed by IWetterichl {2004) as a function of 
the amount of DE at z = 0, the present EoS parameter, w , and 
an average value of the energy density parameter at early times, 
il, A . The EoS parametriza tion that we choose is the one studied 
bv brossi & Springell ( 120091) : 



w(z) 



(1 + C ln(l + z)) 2 



In the above relation the quantity C is given by 
3 wo 



(2) 



(3) 



and characterizes the redshift at which a constant EoS turns into a 
different behaviour according to the presence of DE at early times. 
Since both EDE and ACDM models have to reproduce the observed 
cluster abundance at low redshifts, in the EDE model we expect 
structures to form earlier and to have slower evolution of the halo 
pop ulation that i n the ACDM one. 

AlametaU d201ll ) use d the EDE parametrization by 

ICorasaniti & Copelan j J2003I) to forecast constraints on these 
models from the abu ndance of X-ray c lusters expected in the 
eROSITA survey (e.g. IPredehl et alj|2007h and fr om the SZ power 
spect rum from the South Pole Telescope (SPT. IStaniszewski et al.l 
2009). In our analysis, we use a different parametrization of EDE 
models, we derive forecasts for high-sensitivity X-ray surveys, 
and, particularly, we include the constraints from the redshift space 
power spectrum of clusters. 



a/3 ■ 



d 2 ln£ 



(4) 



where £. is the likelihood of the observable (e.g. Dodelson 2003). 
This can be used to understand how accurately we can estimate 
the value of a vector of parameters p for a given model from one 
or more data sets, under the assumption that all parameters have a 
Gaussian distribution. 

Following th e appr oach of iHolder et al.l (2001) and 
Majumdar & Mohr (2003), the FM for the number of clus- 
ters, Nij„, within the /-th redshift bin and m-th bin in observed mass 
M ob , can be written as 



z 



dN Lm dN^ 1 
dp a dp p N u , 



(5) 



where the sums over / and m run over redshift and mass intervals, 
respectively. In this notation, M? b „ = M thr (z), where M thl (z) is 
defined as the threshold value of the observed mass for a cluster to 
be included in the survey. We write the number of clusters expected 
in a survey with a sky coverage AO., with observed mass between 



<,andM?» +1 , 



and observed redshift between z, and rff. as 



Nl, 



AO. 



r 



dz- 



cIV 



dzdQ. 
"<n dW b 



fv+i dM ob r°° . 

— - dMn(M,z) P (M°"\\M). 
Jm? m Jo 



(6) 



In the equation above, dV/dz is the cosmology-dependent comov- 
ing volume element per unity redshift interval and solid angle, and 
n( M, z) the halo mass fu nction. We assume the expression provided 
bv I Jenkins et all (2001) for the mass function of halos, with mass 
computed at the virial overdensity, A v j r = 324, for the cosmologi- 
cal model assumed in their simulations (see their equation (B4)). 
We verified in test cases that all our results are left unchanged 
if we use instead the more recent c alibration of the mass func- 
tion proposed bv iTinker et all {2008). Furthermore, p(M ob \\M) is 
the probability to assign to each cluster with true mass M an ob- 
se rved mass M ° b . Thi s probability is defined using the prescription 
of lLima&HvJ (2005), which takes into account the presence of 
a lognormal-distributed intrinsic sca tter in the scaling r elation be- 
tween observable and mass (see also lSartoris et alj|2010l . S10 here- 
fater): 



„ h exp[-x 2 (M°*)] 
p(M ob \\M) - 1 



where 



x(M ob ) 



InM" 



(7) 



(8) 



Here B M is the fractional value of the systematic bias in the mass 
estimate and <r lnM is the intrinsic scatter in the relation between the 
true and the observed mass. A negative value for B M corresponds to 
a mass underestimate and, therefore, to a smaller number of clus- 
ters included in a survey, for a fixed selection function. The intrinsic 
scatter has the effect of increasing the number of clusters included 
in the survey. In fact, the number of low-mass clusters that are up- 
scattered above the survey mass limit is always larger than the num- 



© 0000 RAS, MNRAS 000, 000-000 



4 Sartoris et al. 



ber of rarer high-mass clusters wh ich are down- scattered below the 
same mass limit (e.g.. lCunh3l2009l and references therein). 
Including equation 10 into equation {6]l, we obtain 



Nr, 



A£l pft 

I 



dz 



dV 
dzdQ. 



dMn(M,z) [erfc(x m ) - erfc(x m+1 )] 



(9) 



where erfc(x) is the complementary error function. 

In equation ((9}, we assume the sky coverage Afi to be in- 
dependent of the limiting mass threshold or, equivalently, of the 
cluster flux. This formalism can be easily generalized to include a 
flux-dependent sky coverage, which is due to the general degrada- 
tion of resolution and sensitivity over the the field of view of X-ray 
telescopes. Finally, we assume that errors on the cluster redshift 
mea surements can be ignored (see discussion in Section [3~T1 see 
also lLima & Hull2007l for a presentation of a method to include the 
effect of redshift errors in the computation of the Fisher Matrix for 
cluster number counts). 



2.2 Power spectrum 

A new aspect in this analysis, with respect to one presented in S10, 
is the inclusion of the distorted anisotropic mapping between the 
real space density field and measurements in redshift space caused 
by pec uliar velocity (redshift space distortions). Following Kaiser 
( 1987), the redshift space matter power spectrum P (k,fi) in the lin- 
ear regime acquires a dependence on the cosine of the angle be- 
tween the wave number k and the line-of-sight direction, \i, ac- 
cording to 



P(k,fi,z) = (beff + fvff D 2 (z) P (k) 



(10) 



Here, P(k) is the matter power spectrum in real space, b e ff is the lin- 
ear bias weighted by the mass function (see equation (20) in S10), 
and f(a) = dlnD(a)/dlna is the so-called growth function, i.e. the 
logarithmic derivative of the linear growth rate of density pertur- 
bations, D(a), with respect to the expansion factor a. Here we as- 
sume for the bia s of h alos of mass M the expression provided by 
ISheth&TormerJdl999l) . 

The average cluster power spectrum calculated within a given 
redshift interval, P ch can then be written as 



r dz%~n\z)p(n,k,z) 

Pl^(n,k,z0=' 



^dz%n\z) 



(11) 



where h = f Q dMn(M,z) [erfc(x m ) - erfc(x m+1 )]. Therefore, the 
Fisher Matrix for the power spectrum of galaxy clusters is 

1 v-i d\nP c i(iij,k m ,zi) 3 In 
= ft? L Qf a JZ V lmi k m AkAMl2) 



i.m.i 



(e.g. lTegmarkll997l ; lFeidman et alJl994l ; lAbramcl201ll) where the 
sums over /, m, i run over bins in redshift, wavenumber k and an- 
gle fi, respectively. The quantity V cff (pi, k, z) in equation d 12b is the 
effective v olume accessib l e by the survey at re dshift z, at wavenum- 
ber k (e.g. lTegmar k 1997; S artoris et alJuOlOh . 

For the matter power spectru m, we adopt the expres sion for 
the Cold Dark Matter provided bv lEisenstein & Hul dl998l) . which 
includes the effect of BAOs. In order to quantify the information 
carried by the BAOs detection, we use also the power spectrum 
shape from lEisenstein & Hul d 19981) . which smoothly interpolates 



through the oscillations. Moreover, we study the geometric infor- 
mation carried by the shape of the power spectrum by describing 
it with a general free parameter T, and thus ignoring its CDM 
specific relation T = Sl m h. 

In our analysis, we assume the following reference values 
for the cosmologic al parameters, consis tent with the WMAP-7 
best fitting model dKomatsu et al.l 1201 ll) : Q.„, = 0.28 for the 
present-day matter density parameter, Q.^ = for the contribution 
from the curvature, cr 8 = 0.81 for the normalization of the 
power spectrum, £1\, = 0.046 for the baryon density parameter, 
Ho = 70kms~ 1 Mpc~ 1 for the Hubble parameter, n = 0.96 for 
the primordial spectral index. For the DE EoS parametrization of 
equation ([TJ, we take wq = -0.99 and w a = as reference values, 
while for the EDE model of equation I0 we assumed the reference 
values of w = -0.93 and n crfc = 2 x 10~ 4 . Therefore, we have 
in total eight cosmological parameters, which are left free to vary 
in the computation of the number counts and power spectrum 
Fisher Matrices denned in equations l[5} and lQ/2}- We note that for 
the above choice of the cosmological parameters, the reference 
DE model of equation (Q3 is consistent with the WMAP-7 results 
on CMB anisotropies. For the reference EDE model, we adopt 
the same values of the non-DE parameters, including erg. This 
implies that both models are chosen to have the same low-redshift 
normalization as to provide the same cluster number counts 
(see also Fig. Q3, instead of being normalized to CMB. In the 
following, constraints on cosmological parameters will be shown 
in the (fl,„,<x 8 ) and (w ,w„) planes, while marginalizing over the 
remaining parameters. 

Unless otherwise stated, all the constraints that we present 
in the following include the prior information expected from the 
measurement of the CMB anisotropies with the Planck experiment. 
This prior probability has been computed for each of the two ref- 
erence DE models based on equations (T) and {2}. Cosmologi- 
cal constraints from Planck are derived by followi ng the descrip- 
tion presented in the DETF lAlbrecht et all |2009) and by using 
the method described in Ra ssat et"aT1 d2008l) . We conservatively as- 
sume that only the 143 GHz channel will be used as science chan- 
nel. This channel has beam width 9 fwhm = 7.1 arcmin and sensi- 
tivities cr T = 2.2fiK/K and o> = 4.2/iK/K. We take / sky = 0.80 
as the sky fraction in order to reduce the impact of galactic fore- 
grounds. We use the minimum /'-mode, £ min = 30 in order to avoid 
problems with polarization foregrounds. We choose as fiducial pa- 
rameters = (a> m , 9s , In As , a>b, «s , t), where 9s is the angular size 
of the sound horizon at last scattering, In As is the logarithm of 
the primordial amplitude of scalar perturbations and t is the opti- 
cal depth due to reionization. After marginalizing over the optical 
depth, we calculate the Planck CMB Fisher matrix in the param- 
eters (fi,„,fl rf( ,,/j, <T 8 ,n fc , w,«s) by using the appropriate Ja cobian 
of the involved parameter transformation dRassat et al.l l2008). Here 
w is a two-component vector which includes the parameters of the 
two DE models considered here: w = (w ,w a ) for the model of 
equation l[T) and w = (wo, £i e ,de ) for the model of equation {2j. 



3 ANALYSIS AND RESULTS 
3.1 Characteristics of the surveys 

In X-ray flux-limited cluster surveys, when the cluster redshift is 
known, the flux limit for cluster detection can be translated into 
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Table 1. Characteristics of the envisioned WFXT surveys. Sky coverage SI 
(in sq.deg.) (Column 2). Flux limits in the [0.5-2] keV energy band (units 
of 1 0~ 14 erg s -1 cm~ 2 ) for clusters detected with ~ 30 counts (Column 3, see 
also lGiacconi et alj|2009t) and with >1500 counts (Column 4). Number of 
clusters, corresponding to /1500. detected at z > (Column 5) and z > 1 
(Column 6). 



Survey 


n 


/det 


7l500 


Ni50o(z > 0) 


Wl500(z > 1) 


Wide 


20000 


0.5 


15.0 


8471 





Medium 


3000 


0.1 


3.0 


8435 


220 


Deep 


100 


0.01 


0.3 


1509 


375 



a mass limit, based on a previously calibrated relation between 
X-ray luminosity and mass. In order to account for the uncertain 
knowledge of this relation, the so-called self -calibration method 
had been proposed by different auth ors (e.g. iMaiumdar & Morirl 
l2003l ; lLima & Hull2005l : ICunhdl2009l S 10). In this approach, the 
uncertainty in the relation between mass and observable for very 
large samples can be described by an intrinsic scatter and a system- 
atic bias in the estimate of cluster masses. Thus, parameters defin- 
ing the scaling relation are treated as fitting parameters to be deter- 
mined along with the relevant cosmological parameters. Clearly, 
a more realistic a nd direct method to es tima te cluster masses is 
the one adopted in Vikhli nin et all J2009h and lMantz etafl J2010h 
using relatively small samples of clusters based on ROSAT data. 
Deep follow-up Chandra observations with more than 10 3 photons 
for each cluster allowed them to measure mass proxies which are 
closely related to cluster mass, with small (IS 10%) intrinsic scatter. 
Examples of such proxies are the total gas mass, M gas , calculated 
within a fiducial aperture radius, or the product o f gas mass and 
tempe rature, Y x = M glis Tx, originally introduced bv lKravtsov et al.l 
fcOOfj) . 

Next generation of cluster surveys, with high-sensitivity and 
good angular resolution such as WFXT, will yield large subsamples 
of clusters, each detected with a large enough number of photons 
to enable direct measurements of these mass proxies. As additional 
key benefit, these high quality data will allow cluster redshifts to be 
measured from the Fe-K 6.7 keV line in the X-ray spectra, without 
resorting to dem anding optical spectroscopic follow-up campaigns. 
lYu et all d20 1 lh carried out a blind systematic search for K-shell 
and L-shell Fe line complex from Chandra data. Using a sample of 
46 clusters in the Chandra archive, they found that the cluster red- 
shift can generally be measured from X-ray data with a precision 
of Az~ 0.01 when at least 10 3 counts are available. 

In this paper, we derive cosmological forecasts for three ref- 
erence WFXT surveys, which are complementary in terms of sen- 
sitivity and sky coverage: a wide survey covering 20,000 sq.deg. 
down to a flux in the [0.5-2] keV band of /1500 = 1-5 x 
10~' 3 erg s~'cnr 2 , a medium survey covering an area of 3000 
sq.deg. and reaching / 1500 = 3 x 10~ 14 erg s~'cirr 2 and a deep sur- 
vey over 100 sq.deg. with /1500 = 3 x 10 ~ 15 erg s^'cirT 2 (see Ta ble 
QJ. Given the WFXT performances (e.g. iGiacconi et al]|2009l) . at 
these flux limits, clusters are detected with at least 1500 counts, 
thus allowing precise measurements of robust X-ray mass proxies 
and Fe-line based redshifts. 

To convert cluster fluxes to masses, we follow the same pro- 
cedure described in S10 and we refer to that paper for more de- 
tails. We use the relation b etween X-ray luminosity and M500 cali- 
brated bv lMaughanlj2007l) . who analysed Chandra data to estimate 
masses from Y x for an extended sample of clusters over the redshift 
range 0.1 < z < 1.3. 



In the following (see Sect. l3.3.Tt . we show the impact on cos- 
mological constraints of setting strong priors on the parameters 
defining the relation between mass and X-ray observable, when at 
least 1500 photons per cluster are available to measure Y x or M gas 
proxies, as well as the cluster redshift. 



3.2 The mass-observable relation 

Besides the eight cosmological parameters, our Fisher Matrix anal- 
ysis also constrains the parameters which specify the redshift de- 
pendence of the fractional mass bias, B M , and o f the intrinsic scatter 
o~inM- Since curre nt data and simulations (e.g. IStanek et al.ll2010] : 
Fabia n et alj|20 111) show no significant evidence for a mass depen- 
dence of these parameters, in the following we do not consider this 
possibility. According to S10, we assume: 

B M (z) = B m {\+zf 
o-iniu(z) = cr lnM0 (l + zf . (13) 

In this way, we have four parameters, B M0 , o- lnMi0 , a and B, which 
account for the uncertain knowledge in the relation between ob- 
servables and mass (we refer them to hereafter as mass parame- 
ters). We consider a reference value B M0 = -0.15 for the mass bias 
at z = and a = for its evolution. This value of B M .o implies that 
X-ray masses are assumed to be underestimated by 15 per cent. 
This is in line with the level of violation of hydrostatic equilibrium 
found by different autho rs from the analysis of hyd rodynamic sim- 
ulations of clusters (e.g. lBorgani & Kr avtsov 2009, and references 
therein), and from the comparison between weak-lensing and X-ray 
masse s (e.g.. iMahdavi et alj|2003 : IZhang et al]|2010l : lokabe et all 
2010). We also assume an intrinsic scatter <T\ nM $ = 0.25, with 
B = fo r its evolu t ion, c onsistent with the M 500 -L x relation mea- 
sured bv lMaughanl J2007t) . We refer to S10 for a more detailed dis- 
cussion on the choice of these parameters. Following [Lima & Hil 
we point out that we use the variance crf nMQ and not the scat- 
ter as the varying parameter in our Fisher matrix analysis. In fact, 
this variance controls the excess of clusters which are up-scattered 
above the selection threshold, with respect to those that are down- 
scattered. 

In summary, we have four mass parameters that add up to the 
eight cosmological parameters for which we compute the Fisher 
Matrix. In order to quantify the effect of the uncertain knowledge 
of the mass parameters, we set in the following four different levels 
of prior. In order of constraining strength, they can be described as 
follows. 

1. No prior: all the four mass parameters are left free to vary by 
assuming no prior knowledge on their range of variation. 

2. Weak prior: we assume AB M ,o = 0.05, Aa = 1, A<x 2 nM0 = 0.2 
and AB = 1 for the lcr uncertainty with which the four mass param- 
eters are assumed to be known. The above value of AB M0 reflects 
the current uncertainty between different calibrations of violation 
of hydrostatic equilibrium from simulations and from the compar- 
ison of weak-lensing and X-ray masses. The reference value of 
Aa allows for a large variation of the uncertainty with which we 
can calibrate this violation as a function of redshift. The values of 
Act? MQ and AB are rather conservative choices, in view of the large 
number of clusters available from future surveys which will allow 
an accurate estimate of the scatter in mass-observable scaling rela- 
tions. We refer to S10 for a further discussion on the choice of this 
prior for the mass parameters. 

3. Evolution strong prior: in order to emphasize the role played 
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Table 2. Prior on mass parameters assumed in the four cases under study 
(see text) and relative Figure of Merit (FoM) [see equation U4H from the 
combination of three surveys. The analysis is carried out by including the 
Planck prior. 

Reference 



values: 


B M ,o = -0.15 


a = 


<M,0 = (°- 25 ) 2 


/? = 


Cases: 


Strong 


Evolution 


Weak 


No Prior 


AB M ,o 





0.05 


0.05 


1 


Aa 








1 


1 


Acr f„M.O 





0.2 


0.2 


1 


A/J 








1 


1 


FoM 


106 


91 


64 


61 



by the uncertain redshift evolution of the mass parameters, we as- 
sume in this case the uncertainty in B M0 and of M0 to be the same 
as in the weak prior case, while we assume their evolution to be 
known to good precision, so that Aa = and A/? = 0. 

4. Strong prior: in this case we consider the uncertainties in the 
calibration of the mass-observable relation are so small to be ne- 
glected. While this assumption is not realistic in an X-ray clus- 
ter survey providing only detection of clusters, it is plausible for 
a high-sensitivity survey which provides measurements of robust 
mass proxies for all the clusters above the flux limits discussed in 
the previous section. 

We summarize in Table ((2) these different choices of the 
uncertainties in the mass parameters for the different priors that we 
assume. 



30 
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ACDM - 
EDE 
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Redshift 



Figure 1. Comparison between the observed redshift distribution of X-ray 
selected galaxy clusters (histogra m with symbols w ith errorbars) from the 
ROSAT Deep Cluster Survey -3 iRosati et alj|2002l 1 1998h and predictions 
from the reference DE model based on equation (TJ (red solid curve) and 
EDE model of equation (2) (green dotted curve). Errorb ars on observ ational 
data points correspond to lcr Poissonian uncertainties (Gehrels 1986). 



Before proceeding with the derivation of forecasts for con- 
straints on cosmological parameters, we verify that our fiducial cos- 
mological models, with the above reference choice for the mass pa- 
rameters, match available observational data on X-ray cluster sur- 
veys. To this purpose, we show in Fig. [T] a comparison between 
the predicted and the observed redshift distribution for the ROSAT 
Deep Cluster Survey to an X-ray flu x limit of 3 x 1(T 14 erg s^ cm -2 
in the [0.5-2] keV band (RDCS-3, iRosati et aljfl998l . l2002h . We 
stress that this is not meant to be a fit to an observational measure- 
ment of the cluster abundance up to z ~ 1, but rather a test that 
our reference model is consistent with current observations. The 
redshift distributions for the two reference DE models have been 
obtained by convolving the predicted redshift distributions with the 
flux-dependent RDCS sky coverage, which provides complete in- 
formation on the survey selection function. We adopt a minimum 
luminosity of L x [0.5 - 2keV] = 10 42 erg/s which is appropriate for 
the RDCS selection function. 

The good agreement with the data indicates that the our refer- 
ence model can be used to provide a realistic extrapolation of the 
evolution of the cluster mass function over redshift and mass ranges 
which are not probed by currently available data. 

In Fig.|2]and[3] we show the cumulative redshift distributions 
of clusters to be observed in the Wide, the Medium and the Deep 
surveys according to our reference DE and EDE models, respec- 
tively. In the lower panel of Fig. [3] we also show the ratio be- 
tween the cumulative redshift distributions obtained for these two 
DE models by combining the three surveys. In order to reproduce 
the observed abundance of clusters at low redshift, in EDE mod- 
els structures start to form earlier and the halo population follows 
a slower evolution than in the ACDM prescription. We note that 
the combination of the three surveys would provide about 2 x 10 4 



clusters with sufficient photons to allow robust measurements of X- 
ray mass proxies. The total number of clusters is dominated by the 
Wide and the Medium surveys. The Deep survey is expected to pro- 
vide a smaller number of clusters at low redshift, due to the smaller 
survey area. However, the number of distant clusters at z > 1 is 
dominated by the Deep survey, owing to its higher sensitivity (see 
also S 10). 

3.3 Cosmological constraints 

Having defined the reference cosmological model and the charac- 
teristics of the X-ray cluster surveys, we present in this section fore- 
casts on constraints of DE EoS parameters. We show our results in 
terms of constraints on the (Q.,„, <x g ) and the (yv , w a ) plane at the 68 
per cent confidence level, after marginalizing over the other cos- 
mological and mass parameters, and in terms of Figure of Merit 
(FoM). The concept of FoM for DE constraints wa s introduced in 
the D ark Energy Task Force (DETF) report (e.g., lAlbrecht et al.l 
2009) in order to quantify the knowledge on DE EoS parameters 
that future cosmological experiments can reach. In general, the 
FoM for the capability of an experiment to constrain a pair of cos- 
mological parameters (p,, pj) can be defined as 

FoM = (det [Cov(p h /?,)]) " 2 , (14) 

where Cov(pj,pj) is the covariance matrix between the two inter- 
esting parameters. With this definition, the FoM is proportional to 
the inverse of the area encompassed by the ellipse representing the 
68 per cent confidence level for model exclusion. 

In the computation of the cluster number counts Fisher Matrix, 
equation lO, Ni .„ is calculated within intervals of observed redshift, 
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Figure 2. The cumulative redshift distributions of clusters detected with 
more than 1500 net counts in the three surveys, as predicted by the ref- 
erence DE model of equation {T}. Dashed (green), dotted (blue) and dot- 
dashed (cyan) curves are for the Wide, Medium and Deep WFXT surveys, 
respectively, while the solid (red) curve represents the sum of the three. 

with width Az = 0.05 out to z ma x = 2. As for the observed mass, 
we use bins of width AlogM = 0.01, extending from the lowest 
mass limit determined by the survey selection function at a given 
redshift, up to 10 16 /i~'M o . We have verified that with this tight bin- 
ning in mass we saturate information provided by cluster number 
counts to constrain cosmological and mass parameters. 

In the computation of the power spectrum Fisher Matrix, given 
by equation dl2t . the average cluster power spectrum denned by 
equation 1 11 It is calculated by integrating over redshift intervals 
having constant width Az = 0.2. This binning, which is coarser 
than the one adopted for the number counts, was chosen as a 
compromise between the need of extracting the maximum amount 
of information from clustering evolution and the requ est of lim- 
iting the covariance between adjacent z-intervals (e.g.. IStril et al.l 
l20ld) . Indeed, the contribution from different z-bins can be added 
in Fisher Matrix defined by equation l l 1 2t only if they carry statis- 
tically independent information. Using small redshift bins implies 
that the neighbouring bins are significantly correlated. In this case, 
the covariance terms between different redshift intervals should 
be included in the likelihood function entering in the expression 
of equation i 1 2t for the power spectrum FM. As for the wave 
number, we consider a minimum value of k = 0.001 Mpc" 1 ; the 
choice of this minimum value does not change the final results, 
because extremely large wave modes are not sampled by the sur- 
veys used and, therefore, do not provide any contribution to the 
Fisher Matrix. The maximum value chosen is k max = 0.3 Mpc~'. 
This choice derives from the need to maximize the information 
extracted from the three surveys, while avoiding at the same time 
the contribution from small-scale modes where the validity of the 
linear bias model is compromised by the on s et of non-linearity 
(e.g. |Percival& White] |2009| ; IStril et al.1 l2O10l : iRassat et alj|200 8l 
see also S10 for a quantitative analysis of the dependence on 



jTio 2 



z 

"« 1.0 
Q 

A 0.5 

z 





— 1 — 1 — ■ — 1 — ■ — 1 — r 


Wide 












Deep 






Total : 








r 


\ 


■ 




1 ' ' ' ' 1 ' 


» — » — > — t — » — ■ <y 




1 .... 1 .... 1 ... . 



0.0 0.5 1.0 1.5 2.0 

Redshift 

Figure 3. The cumulative redshift distributions of clusters detected with 
more than 1500 net counts in the three surveys, as predicted by the refer- 
ence EDE model of equation (2)- Dashed (green), dotted (blue) and dot- 
dashed (cyan) curves are for the Wide, Medium and Deep WFXT surveys, 
respectively, while the solid (red) curve represents the sum of the three. In 
the bottom panel we show the ratio between the cumulative redshifts distri- 
butions for the combined survey, as predicted by the reference EDE model 
and the reference DE model shown in Fig. [2] The shaded region represents 
the Poissonian error on this ratio. 

k m ax °f FM constraints for n on-Gaussian models). In particular, 
ICrocce & Scoccimarrd d2008l) studied the non-linear evolution of 
BAOs in the Dark Matter power spectrum and correlation function. 
They showed that at k = 0.18 Mpc~' the power spectrum predicted 
by the linear theory is lower by a factor of 1.2 at z = with re- 
spect to the non-linear power spectrum. However, we stress that 
the contribution of information to the Fisher Matrix carried by the 
power spectrum at differ ent redshifts decrea ses for both very high 
and very low values of k dSartoris et alfeoioh . The contribution of 
the power spectrum directly depends on the effective volume, (V e ff 
in equation [T2l. This quantity depends on the power spectrum it- 
self, which is set by the bias parameter, and on the level of Poisson 
noise, which is set by the number density of clusters and is maxi- 
mized at k = O.OlMpc -1 (see Fig.7 of S10). 

Wavenumber bins have been chosen to have log uniform 
width Alogfc = 0.1. Lastly, introducing redshift space distortions 
information, the power spectrum acquires a dependence on /j, 
which is defined as the cosine of the angle that k makes with the 
line of sight (equation U0\)- This implies that the Fisher Matrix 
also involves a sum on /u that runs from ji = -1 to fi = 1. We 
choose to divide the interval of ji into 9 bins. An odd number of 
bins is required so that the central bin samples the values of the 
redshift space power spectrum computed for /z = 0. In fact, this 
number of bins maximizes the power spectrum at different redshifts 
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Table 3. Figure of Merit and r.m.s. uncertainty in the DE EoS parameters, 
for the WFXT three surveys, and for their combination. The analysis is 
carried out by including the Planck prior and assuming strong prior for the 
mass parameters. 



Surveys 


Deep 


Medium 


Wide 


Total 


FoM 


20 


60 


17 


106 




0.20 


0.097 


0.14 


0.064 




0.94 


0.54 


0.70 


0.41 



and wavenumbers and, therefore, the contribution of V e ff (Eq|12t 
to the Fisher Matrix. We also verified that a larger number of 
bins does not tighten constraints from the redshift space distortions. 

The results of our analysis are presented in Fig. [4] where we 
plot the 68 per cent confidence levels on the (fl„,,<x 8 ) and (wo,w a ) 
plane, in the left and right panels, respectively. In each panel, we 
show the contours obtained for each of the three surveys and for 
their combination. Contours are all obtained by combining infor- 
mation from number counts and power spectrum, also including 
the prior information from Planck. A strong prior is also assumed 
for the knowledge of the mass parameters (see Section fOl ). 

The results in Fig. [4] show the trade-off between surveys area 
and depth in constraining different cosmological parameters. As for 
the results on (fl„,,<x 8 ), there is no continuous trend in the con- 
straining power of the three surveys as we reduce the covered area 
and increase sensitivity. The Medium Survey is in fact the one with 
most constraining power, especially for <x 8 , while the Deep and 
the Wide Surveys are somewhat less constraining. Furthermore, the 
three surveys provide comparable constraints on Q. m . This is con- 
sistent with the expectation that constraints on this parameter are 
mainly provided by information on the CMB anisotropies, carried 
by the Planck prior. As for cr 8 , we remind that this parameter deter- 
mines the timing of structure formation, therefore, constraints on 
its value are sensitive to both the number of massive clusters in- 
cluded in a survey, and on the effective redshift range covered by 
the survey itself. In this respect, the Medium Survey provides the 
best compromise between the number of massive clusters detected 
within its area and depth. 

As for the constraints on (wo,w„), we note that their depen- 
dence on the survey area/depth is different from the case of the 
(n m ,crg) parameters. While the Medium survey is still the most 
constraining one, we note that the Deep Survey predict a tighter de- 
generacy between w and w a than the Wide survey. This translates 
into tighter constraints on the redshift evolution of the DE EoS, if 
a prior knowledge on Wo is available, consistent with the fact that 
the Deep Survey covers a larger redshift interval and provides the 
larger amount of clusters for redshift greater then 0.9. This example 
illustrates that the choice of the survey strategy depends in principle 
on the cosmological parameters that one is mostly interested in. In 
Table (O we report the values of the FoM and the r.m.s. uncertainty 
in the DE EoS parameters for each survey and for their combina- 
tion, after marginalizing over the other parameters. The values of 
the FoM in this table confirm that Medium Survey alone carries 
most of the contribution to the FoM obtained by combining the 
three surveys. 

3. 3. 1 Effect of mass parameter priors 

As a first test, we present the effect that using progressively stronger 
priors on the mass parameters has on cosmological constraints. The 



results of this analysis are shown in Fig. [5] where we plot the con- 
straints on cosmological parameters obtained by combining infor- 
mation from cluster number counts and power spectrum, from the 
three surveys together. The Fisher Matrix from the cluster experi- 
ment is also combined with the Planck Fisher Matrix. 

In the left panel of Fig. [5] we show the constraints in the 
(fl,„,cr 8 ) plane. A progressively better knowledge of the relation 
between X-ray observable and cluster mass turns into progressively 
tighter constraints on the cr 8 parameter, while leaving the results on 
n„, basically unchanged. The reason for this behaviour is that con- 
straints on Q.„, are mainly determined by the measurement of the 
CMB anisotropies and by the shape of power spectrum, which how- 
ever only provide rather loose constraints on cr 8 . On the other hand, 
the power spectrum normalization is determined by the growth of 
cosmic structures, which is traced by the evolution of the halo mass 
function. Since a precise measurement of the mass function can 
only be obtained through a detailed knowledge of the mass param- 
eters, it is of little surprise that such parameters determine the ac- 
curacy with which cr 8 can be measured. 

As for the constraints on the DE EoS parameters (see right 
panel of Fig. [5}, we note that improving the knowledge of the mass 
parameters from the no prior (red ellipse) to the weak prior (cyan 
ellipse) case only brings a modest enhancement of the constraining 
power of the surveys. The main reason for this is that constraints 
on DE parameters are here mainly contributed by the evolution of 
linear perturbation growth. On the other hand, constraints on the 
growth are rather degenerate with the uncertainty in the redshift 
evolution of the mass parameters, which is assumed to be rather 
generous also in the weak prior case. Indeed, a more significant im- 
provement in the constraints on the DE parameters is obtained for 
the evolution strong prior case (blue ellipse), which assumes a pre- 
cise knowledge of the parameters determining the evolution of the 
mass-observable relation, with only a slight further improvement 
in the constraints obtained for the strong prior case (green ellipse). 
These results confirm the importance of accurately calibrating the 
evolution of the scatter and bias parameters by measuring different 
mass proxies in high redshift clusters selected in the Deep survey. 
Indeed covering at a high sensitivity an even small sky area allows 
one to obtain a robust calibration of the scaling relations between 
the cluster mass and X-ray mass proxies over a large redshift base- 
line. As shown in Fig. [2] this survey will provide about 400 clusters 
at z > 1 for which measurements of redshift, M gas and Y x will be 
possible, out of which about 100 are expected to lie at z > 1.5. In 
Tableff] by showing the value of the FoM obtained, we summarize 
the results obtained in our analysis on the effect of uncertainties on 
the mass parameters in the determination of DE EoS parameters. 

3.3.2 Combining cluster number counts and power spectrum 

We discuss now how the combination of number counts and power 
spectrum information enhances cosmological constraints. To this 
aim, we show the improvement on constraints obtained by adding 
progressively information from the cluster number counts, the 
mean cluster power spectrum analysis and the CMB prior from 
Planck. The results are presented in Fig. [6] in the (fi„,,cr 8 ) (left 
panel) and (wo, w a ) planes. Constraints are obtained by combining 
information from the three surveys together and assuming strong 
prior on mass parameters. The redshift evolution of the cluster 
number counts sets the direction of degeneracy for the constraints 
on Q.„, and cr 8 . Such constraints are mainly placed on the linear 
growth factor of density perturbations through the mass function. 
Furthermore, since the density parameters contributed by matter 
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Figure 4. Constraints at the 68 per cent confidence level on the (n m ,<xg) parameters (left panel) and on the (wo,w a ) DE parameters (right panel). In each 
panel, forecasts for the Deep, Medium and Wide WFXT cluster surveys are shown with the cyan, blue and green ellipses, respectively, by combining number 
counts and power spectrum information. The red ellipse are the constraints obtained from the combination of the three surveys. All constraints are obtained by 
assuming a strong prior on the knowledge of the mass parameters and combining the Fisher Matrices for cluster number counts, cluster power spectrum and 
CMB Planck experiment. 




Figure 5. Constraints at the 68 per cent confidence level on the (tl m , o-g) (left panel) and on the (wo ,w a ) (right panel) planes. Different ellipses correspond to 
different assumptions on the prior for the mass parameters (see text): no prior case (red), weak prior case (cyan), evolution strong prior case (blue), and strong 
prior (green). All constraints are obtained by combining cluster number counts and power spectrum information for the combination of the three surveys. 
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and DE also affect the expansion history of the Universe, we ex- 
pect their values to be constrained by the cluster number counts, 
through the redshift evolution of the comoving volume element. 
The power spectrum analysis provides information on the growth 
rate of cosmic structure through the bias factor, and the RSDs ef- 
fect. Moreover BAO features, that depend on the expansion history 
of the universe (see Section [333J, and the power spectrum shape 
(see Section l3.3.4t are also sensitive to the underlying DM distri- 
bution. 

In Fig. (|6), we show that adding the power spectrum infor- 
mation to the number counts substantially shrinks the contours. In- 
cluding the information from the Planck prior (red contour) further 
contribute to tighten the contours in the (fl m , erg) plane. In order 
to verify whether CMB add information only by constraining the 
curvature of the Universe, we also show with the green contour the 
effect of assuming instead a flat Universe on the cluster constraints. 
In this case, results on (fi,„,<x 8 ) are not drastically improved with 
respect to the case in which curvature is a free parameter, while they 
are significantly worse than with the Planck prior. The reason for 
this result is that CMB anisotropies provide constraints not only on 
the curvature, but also on the Hubble parameter h, on Sli, and on the 
primordial spectral index n s . All these parameters enter in defining 
the shape of the power spectrum, along with O,,,. Therefore, their 
precise determination from the CMB turns into a significant im- 
provement of constraints on the density parameter from the shape 
of the power spectrum. 

As for the constraints on (w ,w a ), their direction of degen- 
eracy changes as a specific geometry is assumed. Imposing the 
flat prior corresponds to fix the redshift at which DE component, 
fl rfe , starts dominating over O m and, therefore, breaks the degen- 
eracy between wo and w a . By including the Planck prior, instead 
of assuming flatness, has a smaller impact than for the (0„,,cr g ) 
constraints. Therefore, even though CMB alone does not provide 
in itself stringent constraints on the DE EoS, it is quite effective 
in improving the corresponding constraints from cluster number 
counts and power spectrum, owing to its leverage on the geome- 
try of the Universe. The Deep and Medium surveys dominate the 
cluster counts at z~ 0.4 (see Fig.[2jl, thus improving the constraints 
on the growth rate of perturbations in a redshift range where it is 
sensitive to the DE EoS. 



3.3.3 Information from Baryonic Acoustic Oscillations 

We quantify now the geometrical information brought by the pres- 
ence of BAOs features in the matter power spectrum. BAOs ap- 
pear as wiggles superposed on t he power spectrum of t he domi- 
nant Dark Matter component (e.g. lEisenstein & Hu 1998). The os- 
cillation scale is proportional to the inverse of the sound horizon 
at the matter radiation equivalence. BAOs carry at lower redshift 
the same information that baryonic oscillations in the CMB pho- 
ton power spectrum provide at the last-scattering redshift. The po- 
sition of the wiggles is related to the amount of dark matter and 
baryons. As Q„,/j 2 increases, the first peak is shifted to higher k val- 
ues and, moreover, the valleys and peaks become slightly narrower. 
The amplitude of the wiggles also depends on Sit as the oscilla- 
tions grow stronger as the baryon fraction increases. In this analy- 
sis, we study the constraints on the (w , w a ) DE EoS parameters a s 
obtained by using the transfer function bv lEisenstein & HvJ [l998), 
which includes BAOs, and by using instead the transfer function 
that s moothly interpolates thr ough the oscillations (see equation 
(30) of lEisenstein & Hv]|l998l) . In the latter case, the presence of 



baryons manifests itself only by modifying the overall shape of the 
transfer function. 

We carried out the analysis including and excluding BAOs 
in the shape of the matter power spectrum used to compute clus- 
ter number counts and bias. In order to better appreciate the in- 
formation carried by BAOs, in both cases we do not assume any 
prior on cosmological parameters, while we use strong prior on the 
mass parameters. As expected, constraints from Wide survey are 
those that benefits most from the presence of the BAOs. This is 
mainly due to the fact that this survey provides the best sampling 
of the long wavelength modes corresponding to the most prominent 
first oscillation harmonics. The inclusion of the BAOs analysis in- 
creases the FoM by a factor of 2.1 in this case. On the other hand, 
no significant information on BAOs is provided by the Medium and 
the Deep surveys. The FoM from the Medium survey does not in- 
crease, while the FoM provided by Deep survey increases by a fac- 
tor of 1.2. In fact, the Deep survey seems to convey slightly more 
information on the BAOs. This is mainly due to the higher num- 
ber density of clusters in this survey, which reduces the noise when 
sampling the BAOs. 

3.3.4 Information from the power spectrum shape 

To quantify the geometrical information encoded in the matter 
power spectrum shape, we fit the shape parameter T in our anal- 
ysis, regardless of its dependence on Q m , h and Q.b which is spe- 
cific to the type of Dark Matter. We compare the results obtained 
by assuming Dark Matter to be Cold with those obtained for a gen- 
eral form of Dark Matter. In the latter case T is treated as a free 
parameter. In the former, under the CDM assumption, the shape of 
the transfer function is given by Sl m h, that specifies the size of the 
horizon at the equality epoch, and by the baryon density parame- 
ter. For this reason, the power spectrum shape carries information 
on the cosmic expansion history. Relaxing the CDM assumption, 
other characteristic scales could affect the shape of the power spec- 
trum. For instance, if massive neutrinos provide a contribution to 
the DM budget, the power spectrum is expected to be suppressed 
with respect to the pure CDM scale on scales smaller than the 
characteristic neutr ino free streaming scale (e.g., Hannestad 2010; 
iMarulli et alJl201lT) . Already available data on the evolution of the 
cluster mass func tion have been use d to set interesting constraints 
on neutrino mass jMantz et alfeoiol) . 

In Fig. [7] we show the expected 68 per cent confidence ellipse 
on the (Cl m ,crg) plane, by combining cluster number counts and 
power spectrum information for the three surveys together, when 
leaving the shape T as a free parameter (blue dotted ellipse) and 
when using instead its CDM expression (red solid ellipse). In order 
to elucidate the effect of a preadopted CDM power spectrum on 
these constraints, we fix to their reference values the parameters 
that, along with Q m , determine the power spectrum shape, namely 
Q. b , the Hubble parameter h and the primordial spectral index n s . By 
removing the assumption of a CDM spectrum constraints become 
weaker. In addition, since the effect is more pronounced for the 
matter density parameter, the direction of degeneracy changes in 
the sense of a milder dependence of <r 8 on Q m . 

The shape of the power spectrum is better sampled by the 
largest survey area covering the widest range of scale. This trans- 
lates into a stronger constraint on the Q,„ parameters and conse- 
quently on the DE EoS parameters. Thus, if we relax the assump- 
tion of CDM for the shape of the power spectrum, the FoM of the 
Wide survey decreases by 7 per cent. This decrement is less pro- 
nounced in the Medium and Deep surveys, whose FoM decreases 
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Figure 6. Constraints at the 68 per cent confidence level on the (Sl m ,crg) parameters (left panel) and on the (wo,w ) DE EoS parameters (right panel). 
Contours, in order of decreasing area, are obtained by including the Fisher Matrix from cluster number counts only (cyan ellipse), adding cluster power 
spectrum information (blue ellipse), further assuming a flat Universe (green ellipse) and adding priors from the Planck experiment while leaving geometry free 
(red ellipse). All constraints are obtained by combining information from the three surveys and assuming the strong prior on the mass parameters. 



respectively by 4 and 5 per cent, owing to their weaker sensitivity 
to the spectrum shape. 



3.3.5 Information from redshift-space distortions 

In this section, we discuss the effect of including information from 
RSDs in the power spectrum analysis. We remind here that we re- 
strict our analysis to the linear regime, while we do not attempt 
to include the non-linear distortions taking place on small scales. 
In this case, the dependence of the power spectrum on the angle 
between line of sight and wavenumber directions is expressed by 
equation d 10b . The inclusion of the RSDs provides additional infor- 
mation on the linear growth rate of perturbations. This test has been 
amply utiliz ed to constrain cosmic growth from galaxy redshift sur- 
veys (e.g.. lGuzzo et al. IHH; IWang et alj[201ol : iBlake et al.llIoTTl 
and references therein). However, no evidence has been reported so 
far on the detection of such distortions in the clustering analysis of 
galaxy clusters. Also, this information has been never included so 
far in the derivation of forecasts on the constraining power of future 
cluster redshift surveys. 

In Fig. [8] we show constraints on the (h>o, w„) DE EoS pa- 
rameters obtained by either including (blue dotted ellipse) or ex- 
cluding (dot-dashed cyan ellipse) RSDs information in the anal- 
ysis of the cluster power spectrum. Both contours represent con- 
straints derived by combining the power spectrum Fisher Matrix 
from the combination of the three surveys, also including Planck 
prior on cosmological parameters and strong prior on mass pa- 
rameters. Quite remarkably, by including information from RSDs 
DE constraints are tightened significantly, thanks to the additional 
constraints provided on the linear perturbation growth rate. This 
leads to an increase of the FoM by a factor of 2.2. By analysing 
the three surveys separately, if we do not add Planck prior, we find 
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Figure 7. Constraints at the 68 per cent confidence level on the (£2,„, o-g) pa- 
rameters by leaving the shape parameter T as a free fitting parameter (dotted 
blue ellipse) or assuming its CDM dependence on SX m , h and fl;,, (solid red 
ellipse). All constraints are obtained by combining information from the 
three surveys together, including number counts and power spectrum. In 
both cases, we assume here the values of the parameters f2/, and n s to 
be fixed at their reference values (see text). No Planck prior are assumed in 
this case, while strong priors are assumed on the mass parameters. 
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Figure 8. Constraints at the 68 per cent confidence level on the (wo,w a ) 
DE EoS parameters, after including (dotted blue curve) or excluding (dot- 
dashed cyan curve) information from redshift space distortions in the cluster 
power spectrum analysis. The constraints are obtained by combining infor- 
mation for the three surveys together, including the prior information from 
the CMB Planck experiment, and assuming strong prior on the mass pa- 
rameters. 



that including RSDs information enhances the value of the FoM 
by a factor of about 35, 7.75 and 6.8 for the Wide, Medium and 
Deep surveys, respectively. The increasing contribution of RSDs 
with survey depth is due to the fact that tighter constraints are ob- 
tained by extending the redshift baseline over which the evolution 
of perturbation growth is followed. 

We emphasize once again that large surveys of galaxy clus- 
ters do have the potential of conveying cosmological information 
from RSDs. This emphasizes the importance of obtaining precise 
redshift measurements for all clusters included in the survey. 



3.3.6 Constraints on Early DE models 

As a final analysis, we derive now forecasts for the constraints on 
the parameters defining the EoS of Early Dark Energy ( EDE) model 
of equ ation 10, which assumes the parametrization bv lGrossi et al.l 
(2009). In Fig.|9]and[T0] we show constraints obtained on the (wo, 
^e.de) parameters. Cluster number count s and bias are comp uted 
by using the standard mass fun ction by I Jenkins et al.l fcOOlf) . As 
discussed bv lGrossi et al.l {2009), the expression of the mass func- 
tion calibrated on N-body simulation according to ACDM model 
is also a reliable description of the one provided by simulations of 
EDE mo dels, at least as long as DE is homogeneous on small scales 
(see also lFrancis et alj |2009). 

Fig. [9] presents the constraints obtained for each of the three 
surveys and for their combination. They are obtained by combin- 
ing cluster number counts and power spectrum information. We in- 
clude constraints from the Planck prior and assume strong prior on 
mass parameters. The results shown in this figure confirm that the 
Medium survey is the one carrying most of the information on the 



DE EoS. Thus we can extend to EDE models what we found in the 
right panel of Fig.|4]for models of equation Q] where DE influences 
the cosmic evolution at lower redshift. Even if the Deep survey 
should be more sensitive to the EDE thanks to the large number of 
clusters at z > 0.8, a higher FoM is found for the Medium survey. 
This is mainly triggered by the number of total clusters provided 
by the Medium survey even at high redshifts. 

In order to analyse the origin of the constraints on EDE pa- 
rameters, we show in Fig. \10\ how such constraints are progres- 
sively tightened as we add information from cluster power spec- 
trum and Planck experiment to the cluster number counts. The 
latter are expected to provide rather degenerate constraints on 
(wo, ii e .de), which is basically associated to the freedom of choos- 
ing a generic geometry of the Universe. The power spectrum anal- 
ysis brings in addition both information on geometry through the 
shape of the transfer function and extra information on pertur- 
bation growth through RSDs. It is well known that the primary 
CMB, with the exception of the integrated Sachs-Wolfe effect, does 
not provide constraints on the dark ene rgy for a non-flat Universe 
( Albrecht et al. 2006; Bean et al. 2001). However, the Q. m param- 
eter is strongly constrained by the CMB and this in turn tightens 
the cluster constrains on DE, as we point out in Section [3.3.2l (see 
also Figs. [3J. Furthermore, adding constraints expected from the 
Planck experiment causes also EDE constraints to be much im- 
proved, while changing the degeneracy direction. The reason for 
this is that in the EDE scenario, the purely geometrical constraints 
from CMB anisotropics become slightly more important due to a 
non-negligible DE contribution to the total energy density of the 
Universe at z ~ 10 3 . The FoM derived from our constraints in- 
creases by a factor of 27 when we add information from the Planck 
experiment to the clusters analysis. In general, this further high- 
lights that tracing cosmic growth over the widest possible range of 
redshift is requir ed in order to ti ghtly constrain the values of DE 



redsnitt is required m order to tignl 
EoS parameters faa & Viej 2009). 



4 CONCLUSIONS 

In this paper, we presented forecasts on the capability of future 
wide-area high-sensitivity X-ray surveys of galaxy clusters to yield 
constraints on the parameters defining Dark Energy (DE) equation 
of states. We considered the standard equation of state (EoS) pro- 
vided by equation (TJ and the class of Early DE models of equa- 
tion {2^. Our analysis was carried out for future X-ray surveys 
which have enough sensitivity to provide accurate measurements 
of X-ray mass proxies and Fe-line based redshifts for approxi- 
mately 2 x 10 4 clusters, thus extending by more than two orders 
of magnitude the size of the cluste r samples pr esently used to de- 
rive cosmological constraints (e.g., I Allen et al. || 201l|). We used the 
Wide Field X-ray Telescope (WFXT) (e.g., iGiacconi et al.l 120091 : 
iMurrav & WFXT Team! l20ld : iRosati et all l201ll) as a reference 
mission concept along with the Wide (20000 sq.deg.), Medium 
(3000 sq.deg.) and Deep (100 sq.deg.) survey configurations (see 
Table [TJ. We based our analysis on the Fisher Matrix formalism, 
by combining information on the cluster number counts and power 
spectrum, also including the effect of linear redshift-space distor- 
tions (RSDs). This analysis has been carried out with the main pur- 
pose of dissecting the cosmological information provided by geo- 
metrical and growth tests, which are both included in the analysis of 
number counts and clustering of galaxy clusters. The main results 
of this study can be summarized as follows. 

(a) When constraining the parameters of the DE EoS of equation 
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(Q3, we further demonstrate the fundamental importance of having 
a well calibrated X-ray observable-mass relation and, most impor- 
tantly, its redshift evolution. 

We verify that the Figure of Merit (FoM) of the DE EoS increases 
up to 106 when we assume a strong prior on the mass parame- 
ters, as resulting from a precise and robust calibration of the mass- 
observable relation, with respect to the case in which no such prior 
is available (FoM = 61) (see Table l[2j). Such an internal calibration 
can be readily achieved from the same X-ray data by having at least 
~ 10 3 net photon counts for each cluster included in the survey. 

(b) We find that the Medium survey is the one carrying most of 
the constraining power (Table (|3]l), since it is expected to yield 
the largest number of clusters out to redshift z ~ 1. As such, the 
Medium survey shows the tightest constraints on the evolution of 
the DE EoS (<x„ = 0.097 and <r u .„ = 0.54) and the corresponding 
highest Figure of Merit (FoM = 60). The Deep survey, although 
covering a much smaller area than the Wide survey, adds an impor- 
tant contribution to constraining DE parameters (FoM = 20). 

(c) We quantify the increase of the constraining power from the 
three surveys separately and from their combination, by adding pro- 
gressively information from the cluster number counts, the mean 
cluster power spectrum analysis and the CMB prior from Planck 
experiment. We summarize in Fig. [TT] the resulting improvements 
on FoM. The slightly different directions of degeneracy of the con- 
straints in the (w , w a ) parameter space from cluster number counts 
and power spectrum explain why the constraints substantially im- 
prove when we consider the two contribution together rather then 
separately. We verified that adding the CMB information improve 
the corresponding constraints on the DE EoS, mostly as a conse- 
quence of the constraint provided by CMB data on the geometry of 
the Universe (right panel of Fig. [6}. 

(d) We find that RSDs carry important cosmological information 
through the linear growth of perturbations, also in the case of clus- 
ter surveys. Indeed, the DE FoM from the power spectrum analysis 
of the Wide survey increases by a factor 35 when including RSDs, 
while increasing by a factor 7.7 and 6.8 for the Medium and the 
Deep surveys, respectively. 

(e) As for the information carried by the shape of the power spec- 
trum, a smaller increase in the FoM is instead measured when in- 
cluding BAOs. In this case the FoM from the power spectrum anal- 
ysis of the Wide survey increases by a factor of 2, while no signifi- 
cant information on BAOs is provided by the Medium and the Deep 
surveys. Furthermore, relaxing the assumption of CDM and treat- 
ing the shape of the power spectrum as a free parameter reduces the 
FoM by a factor of 1.7 in the analysis of the Wide survey. 

(f) The results obtained for the EDE EoS analysis confirm that 
the Medium survey is the one carrying most of the information on 
the DE EoS. This emphasizes once more the importance of finding 
a good balance in the definition of a survey strategy, between the 
redshift range needed to trace cosmic growth and the survey area. 
By extending the redshift range of the sample and with the abil- 
ity to internally calibrate the observable-mass relation, we expect 
to measure the EDE EoS parameter (il e j e ) with an uncertainty of 
cr Qede = 6.6 x 10- 4 . 



Figure 10. Constraints at the 68 per cent confidence level on Early Dark 
Energy (EDE) EoS parameters including the Fisher Matrix from number 
counts of cluster (dot-dashed cyan curve), adding power spectrum (dotted 
blue curve) and adding prior from Planck experiment (solid red curve). All 
constraints are obtained by combining information for the three surveys 
together and assuming strong prior on mass parameters. 



In order to compare the constraints from the WFXT cluster 
surveys to those expected from other cosmological experiments, 
we also com pare in Fig. | l lithe Fo M expected from WFXT to those 
presented bv lAlbrecht et al.l d2006t) for different large-scale struc- 
ture probes. In the DETF report, they showed that Stage II cluster 
projects (ongoing surveys) provide FoM ~ 2 (Fig. II U when com- 
bined with Planck priors. This analysis was carried out for a generic 
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cluster count survey covering 200 deg 2 up to z„, flv = 2, with the 
simple assumption of a constant mass selection function. Accord- 
ing to Stage IV future experiments, by extending the survey area 
to 20000 deg 2 , the FoM rises in the optimistic configuration up 
to ~ 40 with the contribution of the CMB Planck priors. We point 
out that the DETF analysis did not include the constraints expected 
from the power spectrum analysis. A similar value of the FoM was 
also obtained from an optimistic version of Stage IV project for 
BAO analysis based on galaxy redshift surveys, again including 
Planck priors. The constraining power of these optimistic Stage 
IV experiments is somewhat weaker than that of the WFXT sur- 
veys, with the latter having a FoM larger by a factor of ^ 2. In Fig. 
TTJ we also show for reference the FoM expected for an optimistic 
J Stage IV Weak Lensing experiment, which should reach a value 
of about 300. 

We stress that the above forecasts from the WFXT surveys 
are obtained by considering a subsample of clusters with at least 
1500 net photon counts. With this restriction robust mass and Fe- 
line based redshift measurements can be readily available from the 
same survey data, without resorting to external follow-up calibra- 
tions or observations. As such, the derived constraints should be 
considered as rather conservative since they do not include possi- 
ble information carried by clusters detected with a smaller number 
of photons or any other information to constrain mass from external 
observations (e.g., Sunyaev-Zeldvich fluxes, weak lensing masses 
and optical richness from future surveys). Lowering the flux limits 
of the WFXT surveys by a factor of 30 would still guarantee de- 
tection of clusters as extended sources, without however allowing 
a measurement of redshifts and robust mass proxies. Fig. II llshows 
that by including all the detectable clusters, the FoM increases by 
about one order of magnitude, even by assuming no prior on the 
mass parameters to compensate for the lack of robust mass mea- 
surements. We note that the Wide survey provides the largest con- 
straining power for the DE parameters when we include all clusters 
down to the detection limit. In fact, in this case the Wide survey 
dominates the statistics of clus ters counts out to redshift 1.5 (see 
Figure 3 in lSartoris et aljj2010l) V Clearly, the results obtained from 
all the detected clusters must be considered as optimistic, since they 
rely on the possibility of confirming all these extended sources and 
measuring their redshifts with the aid of large follow-up observa- 
tions. 

In general, our analysis emphasizes that for large cluster sur- 
veys to be really useful for cosmological applications, not only 
large samples are needed but also a robust measurement of mass 
proxies is required for a significant fraction of the cluster sample. 



5 As for the clusters analysis, according to the optimistic configuration, the 
mean of the mass-observable relation and its variance per redshift interval 
of Az = 0. 1 is assumed in the DETF report to be determined up to a level of 
1.6%. 

6 As for the BAOs analysis, the DETF report introduces the <Tf param- 
eter, which describes the scatter in the relation between the true and the 
photometric redshifts, tr 2 = Var(z - z p /,or)/(l + z) 2 - In the the optimistic 
configuration cr/r = 0.01. 

7 As for the WL analysis in the DETF report, the r.m.s. bias a~i n (i +z ) be- 
tween the mean z and photometric redshift for galaxies in ln(l + z) for 
each bin of width 0.15 is assumed to be determined with a precision of 
0.001. Moreover, the shear measurement is assumed to be miscalibrated 
by a factor (1 + f ca f) that varies independently for each redshift bin. It 
is assumed that the calibration factor of each redshift bin has a Gaussian 
prior of width <r(f ca i). In the optimistic scenario this parameter was fixed to 
crtJaa) = 0.001. 



This will be possible with future X-ray surveys only with an ad- 
equate combination of survey area, sensitivity and angular resolu- 
tion. Furthermore, our results also indicates that the optimization 
of the survey strategy depends on the class of cosmological models 
that one wa nts to constrain . For i nstance, in our previous analysis 
presented in Sartoris et al. I d2010h . we showed that the Wide sur- 
vey is best suited to constrain deviations from non-Gaussian initial 
conditions, due to its ability to sample the long wavelength modes 
thus detecting a possible scale-dependence of the bias. This differs 
from the conclusion reached in the analysis presented here, where 
instead we conclude that the Medium survey is best suited to trace 
the growth history of perturbation over a large redshift baseline, 
as required to follow the redshift dependence of the DE EoS. As 
already mentioned in the introduction, models of modified gravity 
and of clustered DE represent an other broad class of models for 
which clusters can provide important constraints ( e.g. jRapetti et al.l 
2009; Sch midt et al]|2009l; lLombriser et ai]|2010h . General predic- 
tions of these models are the scale-dependence of the growth facto r 
of perturbations and of the bias function (e.g., |Parfrev et alj|201 lb . 
Large cluster surveys, such as those considered in this paper, have 
the potential of placing important constraints on such signatures 
of deviations from standard quintessence variants of ACDM, espe- 
cially once a suitable survey strategy is chosen. 

The upcoming results from the Planck mission will much im- 
proved CMB priors based on experimental data. At the same time, 
the eROSITA satellite will provide in few years a wealth of fresh 
X-ray view of the large scale distribution of galaxy clusters and 
the evolution of their population out to z ~ 1. While waiting for a 
first X-ray satellite optimized for surveys, such as WFXT, there is 
no doubt that an extensive follow-up campaign will be required to 
provide a robust mass calibration with independent methods for a 
significant fraction of eROSITA clusters, so as to fully exploit the 
cosmological information contained in such surveys. 
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